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Propagating S=l/2 particles in S=l Haldane gap systems 

Karlo PcncQfl and Hiroyuki Shib J 
Department of Physics, Tokyo Institute of Technology 
Oh- Okayame, 1-1-12, Meguro-ku, Tokyo 152, Japan 
Inspired by the recent experimental on Y2- a; Ca a ;BaNi05, we discuss the dispersion relation of the 
s=l/2 particles in the s=l Heisenberg and VBS model in the limit of small hopping amplitudes. 
The effective s—1/2 edge spins mix with the spin of the impurity resulting in one four-fold and two 
two-fold degenarate bands. We briefly discuss the interaction between the s=l/2 particles arising 
from the background. 
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Recently, the family of the quasi one-dimensional ma- 
terials showing the Haldane gapu has been enlarged by 
the charge transfer insulatom Y2BaNiC>5. In this ma- 
terials the Ni 2+ ions with s=l are surrounded by oxy- 
gens forming an octahedron, NiOg- These octahedra 
are linked and form well separated chains, making an 
ideal one-dimensional s=l antifcrromagnet. Susceptibil- 
ity and neutron scattering measurements have observed 
a spin gap of w 100K (9 meV). The two relevant Ni or- 
bitals are 3d^ z 2_ r 2 and 3d x 2_ y 2. The latter is almost lo- 
calized, whiJe the 3d3 Z 2_ r 2 orbital has finite overlap with 
2p z orbitala of the O. 

The importance of Y2BaNiOs is not only being another 
Haldane gapped system, but by replacing the off-chain 
Y 3+ by Ca 2+ one can dope the chain by holes; thus one 
can introduce carriers into a gapped spin liquid. The 
doped material, Y 2 _ :r Ca ;r BaNi05, has been studied by 
DiTusa et al. (Rcf. |). The addition of carriers lowers 
the resistivity and its temperature dependence cannot 
be described as thermal activation across the charge gap. 
Polarized X-ray absorption experiment showed that upon 
doping, the holes go to the 2p z orbital of the O between 
the Ni 2+ ions. Furthermore, new states with S between 
1 and 3/2 per impurity appeared inside the Haldane gap. 

It is not clear at the moment, what a minimal micro- 
scopic model capable to describe the experiments can be. 
Neither is clear the importance of the disorder caused 
by replacing Y by Ca (the effect of the bond disorder 
was addressed in Refs. ||||). If the disorder can be ne- 
glected, then, as far as the low energy physics is con- 
cerned, the effective model Hamiltonian can be that of 
the s=l Heiesenberg model, where additional s=l/2 ob- 
jects are added (see Fig. |l|). The question of localized 
s=l/2 impurities was discussed in Refs. ||J^. However, 
nothing can exclude that the holes can move, and it raises 
a very interesting question: what are the low energy 
properties of doped Haldane systems? For this reason 
here we discuss a simple Hamiltonian to describe such 
systems. It turns out that even in that simple Hamilto- 
nian, to solve the propagation of one spin s—1/2 object 
is not trivial, and the interaction between the spin 1/2 
objects due to the s=l background is rather complicated. 

A similar problem is the propagation of s=0 holes in 
the Haldane gap systemsa. However, it corresponds to 
a simultaneous hopping of two electrons, which can be 



favorable to one electron hopping only in limited circum- 
stances. 

To describe the interaction between the s—1 spins of 
the Ni 2+ ions, we consider the following rotationally in- 
variant spin Hamiltonian 



Hn = J 



^SjSj+i — [3 (SiSi + iY 



(1) 



We will concentrate on the special cases of (3 = and 
(3 = —1/3. The former is the Heisenberg model for 
s=l, while the latter is the Valence Bond Solid (VBSj) 
Hamiltonian, where the ground state function is knownu. 
Both Hamiltonians have a gap in the excitation spec- 
trum. Since the holes break the valence bonds of the s—1 
host, we get open chains segments between the s—1/2 ob- 
jects (see Fig. |2|). The ground state of an open chain is 
four-fold degenerateE3 in the thermodynamic limit. This 
can be simply explained by noticing that the construc- 
tion of the VBS state of an open chain leaves effective 
1/2 spin objects at the ends of the chain. These edge 
spins can form a triplet and a singlet, and the energy 
splitting for TV-site long chain is ~ exp(— N/£) and it 
disappears for N — > oo, where £ is the correlation length. 
One can actually see these edge,-s±ates, and their exten- 
sion is over several lattice sitestil. For [3 = —1/3, the 
correlation length is In 3, while for the Heisenberg model 
(«7. We denote the triplet wave function by T N ' and 
singlet by Sn, where the superscript ±1,0 denotes the z 
component of the triplet. For N even, the singlet is even 
and triplet is odd in parity, while for N odd the parities 
are reversed. 

Now, let us turn to the Hamiltonian describing the 
effect of s=l/2 impurities. First of all, an impurity at site 
i + 1/2 will change the effective J coupling to J\ between 
the s=l spins at site i and i + 1 in the Hamiltonian in 
Eq. (|l|). Next, there is an interaction between the 1/2 
impurity and s—1 spins (see Refs. H,0): 



Hji — 2J'J^ (Si<x; + i/2 + S,+i<7, + i/ 2 ) 



(2) 



and finally, there are the Hamiltonians describing the ef- 
fective hopping of the holes on the O-sites: 



Hi 



*£4 



i 



Hj 



(3) 



where the factor 2 in front of the J and J' is for con- 
venience and er is the spin operator of the s=l/2 spin. 
The operators P, exchanges the occupation of the site 
i + 1/2 and i — 1/2, i.e. if the hole is on site i + 1/2, 
the Pi will move it to i — 1/2 if that site was empty. 
We do not allow twon s=l/2 objects to occupy the same 
site, as it costs lot of energy. Clearly, the Hamiltonian 
H = Ho + Hj> + Hj + Hj is isotropic. 

In the following we will restrict ourselves to the case 
when the parameters t and J describing the propagation 
of the s=l/2 objects are small compared to the Haldane 
gap, i.e. J of the host s=l system. Furthermore, J\ 
and J' are assumed small as well. In this limit the ener- 
getically large excitation of the s—1 host created during 
the motion can be neglected and we will work in the 
subspace where the wave function of the s=l sequences 
are the singlet and triplet wave function described above. 
Our approach is similar to the variational wave function 
applied in Ref. [7]. 

Given the Hamiltonian and the constraints above, let 
us now turn to the calculation of the dispersion relation 
of a single hole. First, let us construct the trial wave 
function of the L-site periodic chain with one impurity: 



S L ak) = -±=]Te^|5 L a;j + l/2) 



T L( 7fc) = -^^ e ^|T L a;j + l/2), 



(4) 



where \Slo", j+1/2) and |Tj,<t; j+1/2) denote the states 
with hole at site j + 1/2. 

Now it is a good point to say something about the ad- 
dition of angular momenta. The state of the two edge 
1/2 spins and the impurity 1/2 spin has one 5 = 3/2 
and two 5=1/2 representations. Since our Hamiltonian 
is rotationally invariant, we expect one four-fold degen- 
erate band with total spin 5 = 3/2 (quartet) and two 
two-fold degenerate bands with 5=1/2 (dublct). The 
5 = 3/2 representation can be constructed as follows: 



\Qt m ) 



\T+]k) 



\Qt 1/2 ) =y/l/3\T+i k) + v/273 |T°T k) 
\Qk 1/2 ) = VW* \T£ T k) + v/273 \T° L | k) 

\Qk 3/2 ) = \ T Eik), 



(5) 



while the states belonging to the 5=1/2 representations 
are 



\D, 

\D] 



+1/2) 
-1/2, 



'2/3|l£|fc>- 




l/3|T°Tfc) 



(6) 



and 



+l/2\ 
k I 

+ 1/2, 



= \s° L n) 

= \S° L ik) 



(7) 



The parity of the states defined above is different, and 
is determined by the parity of the \Tl) and |5l). For 
example, if we define the reflection operator R so that it 
changes the spin at site j to — j, then R\Q k ) — r\Q k ), 
R\Dk) = r\D k ) and R\D k ) = -r\D k ), where r = -1 for 
L even and r = +1 for L odd. 

The rotational invariance of the Hamiltonian implies 
that the Hamiltonian matrix among \Qk) is diagonal, 
and there are matrix elements between the |Z?fc) and 
|-Dfc) states only. Furthermore, since the parity of the 

is different, the matrix ele- 



states I D 



±1/2, 



and \D 



±1/2, 



ments with some hermitian operator A commutable with 
R will satisfy (D k \A\D k ) = -{D- k \A\D^ k ). These argu- 
ments give 



H\D k ) = (e T + b k )\D k ) + ic k \D k ) 
H\D k ) = -ic k \D k ) + (e s + d k )\D k ) 
H\Q k ) = (e T + a k )\Q k ) . 



(8) 



Here £5 and Et are the energies of the \Sl) and \Ti), 
respectively, and for small values of the interaction pa- 
rameters compared to J they depend only on J' apart 
from the finite size effects mentioned earlier. Then, the 
dispersion relation is 



4(fc) = b k +e T +d k +e s ± 
e Q (k) = s T + a k . 



'b k +s T -d k -e s 



(9) 



We can give more precise statements about the pa- 
rameters using the microscopic model of the hopping, 
Eq. (||). First, let us consider the hopping due to Hf. 
Since Hf is a scalar operator in the spin space, it will have 
purely diagonal matrix elements: (T m ak\Hi \T m a'k) = 
2S mm '6 aa 'thf cos k and (Sak\Hj\Sa'k) = 2S aa 'th® cos k. 



However, Hj is a vector operator in the Hilbert 
space spanned by the s=l spins, and we can use 
the Wigner-Eckart theorem to get the following 
identities: (T m ak\Hj\T m ak) = -imaJlij cos k and 
(T+ 1 i:k\Hj\T° t;fc) = 2v^J/ijCOsfc, furthermore 
(Sa;k\Hj\Sa';k) = 0. Similarly, (T+ 1 1; k\Hj\S T; k) = 
-i2V2Jh I ?smk and (T° ^;k\Hj\S V,k) = i2Jh^ sxak. 
The same arguments can be applied for the matrix ele- 
ments of Hj/ as well. Putting all together, we get 

a k = 2(hft + h Q jJ) cos k + J'g 
b k = 2(h?t - 2h Q jJ) cos k - 2J'g 
Cfe = 2V3/ij Jsinfc 

d k = 2h°tcosk (10) 

where the coefficients hf ,h® ,h® , h j and g depends on 
the wave functions |5l) and \Ti) and their size depen- 
dence is governed by the £. 
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Here we will calculate these coefficients for the VBS 
model and the Heisenberg model. 

VBS model: For the VBS model the task is essentially 
simplified due to the knowledge of how to construct the 
ground state wave function. Actually, for our purposes, 
the following decomposition of the four lowest lying wave 
functions turns out to be useful: 



\T° N 
\T N 
\Sn 



= a N \T+^0) - a JV |T°_ 1 +) - (3 N \S N ^+) 
= "aH^-i-) _ a N\T N _ 1 +) - (3 N \S N -iO) 
, = oln\T^_-t) - c(n\T iv _ 1 0) - Pn\Sn-i-) 
) = v/173 (|3^_ x -) - m^O) + |T"_ 1+ )) 



and similarly 



\1%) = 
\T° N ) = 

\T N ) = 



-a N \QT£_ 1 ) + a N \+'I% r _ 1 }-p N \+SN-i) 
-a Ar |-r+_ 1 ) + a N \+T N _ 1 ) - p N \0S N ^x) 
~a N \-T^_ 1 ) + a N \0T N _ 1 ) - Pn\-Sn-i) 
\S N ) = VT/3 (|-T+_ X ) - |0Tg._ a > + | +T^_ X » , 

where the coefficients a at and [3n are given in Tab. |. ft 
is easy to get this decomposition by inspection. We be- 
lieve that it is trivial to get this result from the transfer- 
matrixE3 representation of the VBS wave function. The 
coefficients exhibit the following remarkable properties: 
2ajy + /3]y = 1, and it means that in this decomposition is 
complete in the Hilbert space spanned by the four VBS 
wave functions. Furthermore, = (2 + 3a^_j) _1 ' 2 , 
which allows us to calculate the coefficients recursively. 

With the help of the wave function decomposition pre- 
sented above, it is straightforward to get the parameters 
for the hopping matrix elements: 



2aj 



h? = Pi 
h Q ,=Pl 



hf = i 

hf = 2a L /V3 



fJ 



2(4 



Pi) 



(11) 



furthermore £5 = — 2Ji/3 and Et = (34 
where the energy is measured from £t(Ji = 0) 



80a| 



Jx/9, 
es(Ji = 



0) . Let us comment here that = 1 means that no walls 
were created during the motion for that special process. 

For infinitely large system and J\ = the dispersion 
relation Eq. (||) is simplified to 



e±(fc) = -tcosk - -J' ± - J 3 J 2 sin 2 k + (t cos k + J') 2 



2 - 4 
£ Q( k ) = -gicosfc + - J' 



(12) 



We show some examples of the dispersion relation in 
Fig. |[ An interesting feature of the dispersion relation is 
that for large values of J the minimum moves away from 
k = n (when t > 0) or k = (if t < 0). It means that 



the holes will be described by a two-band model, which 
can have interesting features. 

Heisenberg model: We also calculated the hopping ma- 
trix elements in Eq. ( |To| ) for small clusters of up to 15 sites 
for the more realistic Heisenberg model, where f3 = 0. In 
that case the correlation length is comparable to the clus- 
ter size and the size dependence of the matrix elements 
is large. We have plotted the different matrix element 
on Fig. (^). Although the size is not large enough to 
get good values for L — > 00 limit, we can conclude that 
hf = -G.28±0.01, h Q j = 0±0.005, hf = 0.81±0.01 and 
hj = 0.55 ± 0.01. Furthermore, from Ref. |^ we know 
that g — a and ss — £t ~ ot 2 J\ for J\ <C J, where 
a = 1.0640. We find these matrix elements to be ~ 20 % 
less then those of the VBS Hamiltonian. 

A few words about the validity of the approach pre- 
sented above. During the motion the spin 1/2 object can 
destroy the hidden AF ordered by creating walls. Taking 
this into account, it would give us corrections of the order 
t 2 1 J and J 2 1 J to the dispersion relation. Also with the 
increasing amplitude of the hopping, the upper bands 
will merge with the continuum of the states above the 
gap. On the other hand, the numerical calculation on 
small clusters shows that the qualitative features of the 
lowest band remains even when the hopping amplitudes 
are comparable with the magnitude of the Haldane gap. 

Now, let us turn to the question of what happens if 
there are more than one s=l/2 impurities? We can follow 
the idea that for small values of the t and J, the states 
above the gap are not excited and it is enough to consider 
the four low-lying states of the open chain for the wave- 
function of the s=l sequence in the wave functions. For 
example, a typical state is 



(13) 



We can think of this wave function as an variational 
ansatz. Then, the interaction between the s=l/2 ob- 
jects in this wave function comes from: (i) the hopping 
matrix elements, which depends strongly on the size of 
the open chains, i.e. on the distance of the nearest holes; 
(ii) the energy splitting of the singlet and triplet states of 
the finite chains between the holes. The interaction due 
to (i) is proportional to the hopping amplitudes t and J 
itself, while the strength of (ii) is determined by the J 
of the host system. These effects depend very much on 
the correlation length of the s=l system, and they are 
the smallest for the VBS model. Actually in that case 
the energy splitting is zero and the interaction is due to 
(i) only. These interactions can in principle result in a 
bound state, unless the kinetic energy is large enough. 

Despite of the strong constraints involved in construc- 
tion of the variational wave function Eq. (|l3|), it has still 
a substantial freedom and the properties of the system 
with more impurities remains to be solved. 

As far as the experiments done on Y2_ a; Ca 2 ;BaNi05 
are concerned, we have shown that there are states with 
S larger than the s=l/2 of the impurity in the Hal- 



3 



dane gap. Unfortunately, the parameter range for the 
Y2- 2; Ca 2 ;BaNi05 based on simple electronic model seems 
to indicate that the hopping amplitudes are comparable 
with the interaction between the s=l spins, where our 
approach is valid only qualitatively. 



FIG. 1. Schematic diagram of the Ni-O chain with typical 
low-energy level occupation. The hole is on the second O. 

FIG. 2. Part of the s=l chain with two s—1/2 spins. 
Between them the three s=l spins form an open chain. 
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FIG. 3. Dispersion relation for different values of J/t and 
J' — Ji = 0. The dashed line is £q(/c), the lower solid line is 
for e]y(k) and the upper for ej(fc). Energy is measured from 
es = £t- 

FIG. 4. The matrix elements h for the Heisenberg 
model: hf (triangles), (diamonds), hfj (squares) and hf 
(hexagons) from the top to the bottom. 



TABLE I. The coefficients of the VBS wave function de- 
composition 
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